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Summary 

Analytical  stress  solutions  are  presented  for  a  metallic  rocket  motor  case 
reinforced  with  a  prestrained  fibre  overwind,  with  viscoelastic  properties, 
under  both  constant  and  varying  temperature  histories.  The  treatment  is 
In  terms  of  the  "reduced  time”  proposed  by  Schwarzl  and  Staverman,  with 
the  Williams/Landel/Ferry  form  of  the  "shift  function”. 

The  analysis  is  reduced  to  the  solution  of  a  Volterra  integral  equation. 

For  a  constant  temperature  history  and  using  the  "standard  linear  solid” 
representation  of  the  material  properties  of  the  overwind,  this  equation 
can  he  solved  in  closed  form.  For  varying  temperature  histories  the 
equation  is  solved  using  a  finite-difference  technioue  proposed  by  Lee 
and  Roger  and  improved  by  Margetson.  The  effects  of  fibre  relaxation  on 
the  internal  pressure  required  for  the  initiation  of  yield  in  the  motor 
case  are  evaluated  and  discussed  for  a  number  of  temperature  histories 
and  end  conditions. 
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1  INTRODUCTION 

In  an  effort  to  construct  rocket-motor  cases  with  greater  specific  strengths, 
several  attempts  have  been  made  to  utilise  the  higher  strength  maraging 
steels-  This  has  resulted  in  serious  quality-assurance  problems,  often 
associated  with  reduced  ductility,  and,  in  some  instances,  rocket  motors 
have  been  observed  to  fail  in  a  brittle  manner. 

With  the  choice  of  suitable  constituent  materials,  however,  a  composite 
motor  case  can  be  constructed  with  a  specific  strength  which  exceeds  that 
of  any  homogeneous  steel  design.  A  lightweight  alloy  case  overwound  with 
a  high  strength  polyamide  fibre  is  a  particular  example.  Moreover,  if  the 
fibre  is  wound  under  tension  a  compressive  stress  is  induced  in  the  case, 
reducing  the  net  tensile  stresses  developed  in  operation,  and  effectively 
increasing  the  specific  strength.  However,  this  strength-enhancement 
technique  can  create  problems  associated,  for  instance,  with  the  dependence 
of  the  mechanical  properties  of  the  overwind  on  environmental  conditions 
and  sustained  loading^.  Nevertheless,  the  considerable  possible  gains  in 
specific  strength  often  outweigh  these  potential  disadvantages. 

A  recent  extensive  study  of  the  behaviour  of  the  overwound  rocket  motor 
case2-1*  has  included  analyses  of  the  initial  winding  process,  the  effect 
of  changes  in  the  thermal  environment  and  the  behaviour  of  the  motor  case 
on  pressurisatlon.  The  work  has  been  limited  to  linear  elastic  conditions; 
the  effects  of  fibre  relaxation,  which  is  known  to  occur,  on  the  initial 
stress  state  and  subsequent  yield  behaviour  of  the  motor  have  not  been 
considered. 

In  this  report  procedures  are  developed  for  the  stress  analysis  of  a  metallic 
motor  case  overwound  with  a  thermorheologically  simple,  linear  viscoelastic 
fibre.  The  Young's  modulus  of  the  fibre  is  assumed  to  be  time-dependent; 
the  Poisson's  ratio  and  thermal  expansion  coefficients  are  assumed  to  be 
time-independent.  These  assumptions  are  valid  given  the  temperature 
ranges  considered  in  this  work.  The  case  material  is  assumed  to  be 
linear  elastic  up  to  yield. 

The  development  of  the  viscoelastic  analysis  is  presented  and  techniques 
for  solving  the  resulting  integral  equation  are  discussed.  The  manner  in 
which  the  storage  temperature  affects  the  relaxation  rate  is  considered. 
Methods  for  predicting  the  yield  pressure  behaviour  for  various  winding 
configurations  and  relaxation  rates  are  presented,  and  typical  numerical 
results  are  discussed. 

2  CONSTITUTIVE  EQUATIONS 


2. 1  Temperature-Independent  Analysis 

As  a  precursor  to  the  temperature-dependent  analysis,  in  which  variations 
of  mechanical  properties  with  temperature  are  considered,  the  constitutive 
equations  for  temperature-independent  behaviour  are  briefly  discussed. 

The  viscoelastic  constitutive  law  can  be  expressed  either  by  the  differential 
operator  method5  or,  alternatively,  by  the  more  general  representation  of  the 


3 

UN  UNITED 


UNLIMITED 


hereditary  or  Duhamel  Integral6,  In  which  the  strain  Is  treated  as 
consisting  of  successive  increments,  applied  over  a  sequence  of  tines 
0  t'  t,  giving  the  stress  as: 


t 

oft)  -  t(t)E(0)  -  /  jEft-tO  E(tOdt- 
o  dt' 


(1) 


In  this  equation  oft)  and  eft)  denote  the  components  of  stress  and  strain 
respectively  at  time  t,  and  Eft)  Is  the  relaxation  modulus. 

The  form  of  variation  of  the  relaxation  modulus  Eft)  with  time  is  shown  by  the 
central  curve  In  Figure  1;  particular  importance  attaches  to  the  short-term 
value  EfO)  fie  the  “instantaneous”  modulus)  and  the  long-term  asymptotic 
value  F(«)  fie  the  "equilibrium"  or  “rubbery"  modulus). 


2.2  Temperature-Dependent  Analysis 
2.2.1  Constant  Uniform  Temperature  Change 


In  the  preceding  section  no  allowance  was  made  for  changes  in  the  mechanical 
behaviour  of  the  linear  viscoelastic  material  arising  from  possible  differences 
in  the  winding  temperature  and  the  subsequent  temperature.  For  most  suitable 
overwinding  materials  such  effects  are  second  order  and  may  be  neglected7, 
but  there  may  be  applications  in  which  the  mechanical  properties  are 
significantly  influenced  by  temperature  changes  and  these  effects  cannot 
be  neglected.  This  section  deals  with  such  applications  and  it  is  assumed 
that  the  structure  experiences  a  constant  spatially  uniform  temperature  change 
subsequent  to  winding. 

Early  work  by  Leaderman8  revealed  that,  by  adjusting  the  time  scale,  it 
was  possible  to  represent  the  temperature-dependence  of  several  mechanical 
properties  of  a  wide  range  of  viscoelastic  materials  by  means  of  a  series  of 
master  curves.  Such  materials  were  later  to  be  defined  as  "thermorheologically 
simple"9. 

The  existence  of  this  simple  translational  shift  led  Schwarsl  and  Staverman9 
to  propose  the  following  equivalence  relationship  between  the  relaxation 
modulus  and  temperature: 


FfT.t)  -  E(T0,0 


f2) 


where  c  “  t/ax  fT) 
o 


f3) 


In  the  above  equations  T  denotes  the  temperature,  T0  a  reference  temperature, 
t  the  time  elapsed  from  the  completion  of  winding,  and  e  the  "reduced"  or 
"shifted”  time,  rela'.ed  to  the  real  time  via  the  "shift  function"  aT  (T). 
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Since  the  relaxation  curve  (see  Fig  1)  shifts  to  the  left  for  a  temperature 
increase  and  to  the  right  for  a  temperature  decrease,  it  follows  that  for  T 
greater  than  T0,  a^  (T)  is  less  than  unity  and  vice  versa.  It  is  noteworthy 
o 

that  both  the  instantaneous  and  equilibrium  moduli  are  temperature-independent. 

The  relaxation  modulus  at  an  arbitrary  temperature  can  thus  be  expressed  in 
terms  of  the  reference  temperature  behaviour  using  a  new  time  scale  which  is 
dependent  upon  the  particular  temperature.  Expressing  the  stresses  and  strains 
as  functions  of  reduced  time,  the  viscoelastic  constitutive  relationship, 
Equation  1  becomes 


o(  C)  ”  e(OE(0)  "  /  c')d  C 

o 


(4) 


2.2.2  Time-Dependent  Uniform  Temperature  Change 

Since  the  uniform  temperature  considered  in  Section  2.2.1  was  constant,  a 
simple  expression  for  calculating  the  reduced  time  resulted.  It  is  now 
necessary  to  extend  this  analysis  to  deal  with  time-dependent  temperature 
histories.  Again  spatial  temperature  variations  are  assumed  to  he  zero;  the 
size  and  thermal  dlffusivity  of  the  components  under  consideration  arq  such 
that  this  assumption  is  not  unrealistic.  It  is  also  necessary  to  assume  that 
the  shift  function  aT0(T)  at  a  particular  temperature  is  Independent  of  the 
prior  temperature  history.  Experimental  evidence  by  Williams  et  al  10 
suggests  that  such  an  assumption  is  valid. 

With  these  assumptions  it  follows  from  Equation  3  that 


dc(t)  -  dt/aT  (T(t ) ) 
o 


(5) 


Hence  integrating  and  noting  that  at  t  ■  0,  £(t)  “  0 


t 

C(t)  -  /  dt'/aT  CT(t')) 
o  o 


(6) 


In  establishing  a  general  form  of  the  shift  function,  Williams  et  al  10 
found  that  the  empirical  relationship 


l°g10taTo(T>I 


-k,(T-T0) 

k2+T-T0 


(7) 


in  which  k}  and  V-2  are  material  constants,  gave  good  results  for  a  large 
range  of  viscoelastic  materials. 
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3  VISCOELASTIC  ANALYSIS  OF  TUE  REINFORCED  MOTOR  CASE 

For  a  body  subjected  to  a  series  of  external  loads  and  displacements,  the 
stresses  and  strains  at  any  point  must  satisfy  the  basic  requirements  of 
equilibrium  and  compatibility  and  must  conform  with  the  appropriate 
constitutive  relationships.  Moreover,  the  stresses  and  displacements  at 
the  boundary  must  comply  with  the  applied  external  forces  and  displacements. 


The  gradual  decease  in  the  relaxation  modulus  of  the  fibre,  Eg(c),  will 
result  in  a  time-dependent  expansion  of  the  pre-stressed  motor  case.  If 
u(?)  denotes  this  expansion,  then  from  elementary  thin-walled  cylinder 
theory,  the  Initial  strains  in  the  fibres  and  case  are  modified  according 
to  the  relations 


e 


(i,n) 

6f 


(?) 


<i,n) 

Eef 


(0)  + 


1-1,2, 


>n, 


(P) 


(n) 

e0c 


(?) 


(?) 


(i,n)  (n) 

In  these  equations  Egf  (0)  and  egc  (0)  are  the  Initial  strain 
components  for  the  it”  of  n  fibres  and  the  case  respectively,  as  given  in 
Reference  2,  and  5  is  the  reduced  time  defined  by  Equation  6.  The  quantity 
R  denotes  the  mean  radius  of  the  motor  case. 


When  the  composite  motor  Is  subjected  to  a  spatially  uniform  temperature 
change,  thermal  strains  are  induced  in  the  fibres  and  case,  further 
modifying  the  initial  strains  produced  during  winding.  In  general  the 
coefficient  of  expansion  of  the  fibres  (of)  will  be  less  than  that  of 
the  case  (a,.),  and  hence  an  increase  in  temperature  will  increase  the 
compressive  stress  in  the  case  and  vice  versa. 


From  a  previous  thermal  analysis3,  it  can  be  shown  that  the  strains  in 
the  ith  fibre  layer  and  the  case  resulting  from  a  uniform  temperature  change 
AT((;)  at  reduced  time  x,  are  given  respectively  by 

(i,n)  Sf(t) 

e0fAT  + — j — ’  1  -l,2,...,n, 

(10) 

and 


e 


(n) 
0c  AT 


(?) 


acATU)  + 


Me) 

p 


(ID 


These  strains  consist  of  two  components,  (i)  a  component  due  to  unconstrained 
thermal  expansion  and  (ii)  an  interference  component  (df(c)/R  for  the 
fibres,  6C(?)/R  for  the  case)  resulting  from  the  interaction  of  the  two 
materials. 
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(i,n) 

It  follows  that  the  total  strain  eafTOT^O  for  the  fibres  may  he  obtained 
by  adding  Equations  8  and  IP  to  give 


£8f  TOT  (C) 


(l.n) 


*  A*--  /  y\ 


where 


(1?) 


Aef  (?) 


Mil 

R 


(13) 


Similarly  for  the  case,  the  total  strain  is  obtained  by  adding  Equations 
9  and  11,  giving 


(n)  ,  a 

£8cT0T(  5) 


(n) 

£ec  (°) 


+  At^  ?)  +  Oc  AX(  ?) 


where 


(14) 


As  <  C) 
c 


Mil 

R 


(15) 


As  the  fibres  and  case  are  assumed  to  remain  in  contact,  compatibility  of 


radial 

equal. 


(i,n)  (n) 

displacements  requires  that  tgf  g  (?)  and  £ecAT(c)  are 
It  therefore  follows  from  Equations  10,  11,  13  and  15  that 


Aec(?)  “  Acf(t)  +  (<jf  -  oc)AT{?) 


(16) 


Using  Equation  16  the  quantity  ACc(?)  appearing  in  Equation  14  can  be 
eliminated  in  favour  of  Acf(?)  to  give 


(n)  (n) 

E0cTOt(O  “  E0c  (°)  +  6Ef(?)  +  Ofsno 


(17) 


In  Equation  A  the  strain  e( ?)  is  the  mechanical  strain,  ie  the  total 
strain  minus  the  free  thermal  strain.  Making  the  appropriate  substitution 
for  the  mechanical  strain  from  Equation  12  into  Equation  4,  it  follows  that 


(i,n)  .  (i,n) 

°0f  <I>  ■'(<”)  (0)  +  Aef  ( c) ) 

-  «•”><„> 
o  * 

+  Aef  (?'))d?',l»l,2,...,n, 
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(i,n) 

where  om  (t)  is  the  stress  in  the  itn  fibre  layer  at  reduced  time  £. 
Similarly,  the  circumferential  stress  in  the  metal  case  at  time 
(n) 

Z,  °6c  (?)»  follows  from  Equation  1/  and  the  linear  elastic  stress-strain 

relations  in  the  form 


<»),  \ 
09c  (?) 


(n) 

(tec  (o)+Atf(c) 


(  nc  —  fiL.  1  NT  (  T 1 


(19) 


Kadlal  stresses  are  small  and  have  been  omitted  from  these  equations. 


In  the  absence  of  an  internal  pressure,  it  is  an  equilibrium  reauirement 
that  the  sum  of  the  forces  associated  with  these  stresses  (ie  Eouations 
18  and  19)  is  zero3,  ie 


"  (i,n)  (n) 

I  oef  (C)tf  +  oqc  ( ;)tc  -  0 

i-1  (20) 

where  tc  and  tf  denote  respectively  the  case  thickness  and  the  effective 
thickness  of  a  single  fibre  layer.  Substituting  into  Equation  20  from 
Equations  18  and  19  for  the  fibre  and  case  stresses  results  in  the 
following  integral  equation  for  the  unknown  quantity  Atf(c) 


where 


G(0)  Aef  (  c)  -  /  £G(_£jQ ACfU'Wc'  -  f(?) 
o  ?' 


f(?)  ■  E(tc[e^n)(0)^^  -! 


-(af-ac)  AT(  ?) 


(21) 


and  G(?)  -  nEf(c)tf  +  Fctc. 


(22) 


Equation  21  is  a  Volterra  integral  equation  of  the  second  kind5  which,  for 
a  specified  f(^)  and  G  ( <;) ,  has  a  unique  solution  for  Atf(c). 

With  this  solution,  the  change  in  the  mechanical  strain  in  the  case  At^-ft) 
can  be  computed  from  Equation  16.  The  corresponding  stresses  then  follow 
from  Equations  18  and  19.  Details  of  the  solution  procedure  are  given  in 
the  following  section. 


4  SOLUTION  PROCEDURE 

Solutions  of  equations  similar  in  form  to  Equation  21  may  be  obtained 
using  Laplace  transforms.  This  solution  procedure  eliminates  the  time- 
dependence  of  the  system  equations  in  favour  of  the  transform  parameter, 
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thereby  reducing  the  viscoelastic  analysis  to  an  associated  elastic 
problem.  Inversion  of  the  resulting  equations  yields  the  solution  to  the 
original  viscoelastic  problem. 

For  elementary  models,  eg  Maxwell  and  Kelvin  idealisations7,  analytical 
inversion  is  often  possible,  but  these  idealised  solutions  are  usually  of 
limited  practical  interest.  For  many  practical  problems,  the  relaxation 
modulus  can  be  represented  by  the  Dirichlet-Prony  series  11 


Ff(t) 


Ef(-)  +  l  EfJ 
j-1 


(23) 


where  pc(«>)  and  Fj,  denote  the  long  term  modulus  and  the  of  m  Prony 
coefficients  respectively  and  tj  is  the  relaxation  time  associated  with 
the  jt31  Prony  coefficient. 


Appropriate  numerical  techniques  (eg  those  proposed  by  Bellman  et  al  12) 
can  then  be  used  for  inverting  the  Laplace  transforms.  Using  this  approach 
for  the  present  problem,  spurious  oscillations  in  the  solution  were 
observed,  particularly  at  times  corresponding  to  large  relaxation  rates, 
which  could  not  be  eliminated. 


An  alternative  solution  technique  was  therefore  developed  in  which  the 
integral  appearing  in  Equation  21  was  expressed  as  a  sequence  of  analytical 
expressions,  as  originally  proposed  by  Lee  and  Rogers  13  and  subsequently 
improved  by  Margetson  3I*.  If  the  time  interval  0  to  q  is  divided  into  a 
sequence  of  time  intervals  by  the  time  values,  q^,  i  “  1,2,  ...,  L, 
such  that  ci,  ”  0  and  qj_  ”  C»  and  the  fib*6  strain  change  AEf(q)  is  the 
mean  value  taken  over  the  interval,  Margetson  has  shown  that  the  strain 
change  A Ef(c),  for  q  -  q^,  can  be  obtained  from  the  equation 

“  ,  “UL-tL-l^Tj  + 

dtf  (?L>  *  +  l  ntfEf  j  !sj,L-le 

j-1  (24) 


-<CL-CL-lVTj  , 

*iAef(  qL_i)(l-e  )  }/ 

®  -(Cl.-«L-l>/'Tj^ 

{G(0)  -  \  l  ntfFf  jd-e  } 

2j-l 

In  the  above  equation  the  quantity  Sj(l  denotes  the  recurrence  vector 
_(?L-CL-l)/Tj 


sj,L  ‘  Sj,L-le 

+  ?tA£f(V  +  vwi 


(25) 


(1- 


-Ul-?L-i)/tj 
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Once  Aef(ci_)  has  ^een  evaluated,  the  stress  states  for  the  fibres  and 
case  can  be  evaluated.  It  can  be  shown  14  that  the  tine-dependent  fibre 
stress,  Eouation  18,  is  given  by 

°0f(1,n)^L)  -  Kf(0)Aef(Cl)  +  C0f(1,nVo)Ef(cL) 


"  Ij  EfjSJ.L 


i  =  1,2,. ..,n, 


(26) 


where  S<  ^  is  the  recurrence  vector  defined  by  Equation  25.  For  the 
case,  the  corresponding  stress  component  then  follows  from  the  equilibrium 
condition  given  by  Equation  20. 

The  finite  difference  scheme  outlined  above  has  one  notable  advantage 
over  the  original  technique  proposed  by  Lee  and  Roger13,  in  that  the 
use  of  the  recurrence  vector  dispenses  with  the  storage  of  all  strain 
values  other  than  those  obtained  at  the  last  two  time  intervals,  thus 
reducing  the  necessary  computer  storage  and  improving  numerical  efficiency. 
It  is  noteworthy  that  the  finite-difference  scheme  introduces  no 
complications  when  the  time  intervals  are  unequal.  This  enables  larger 
time  Intervals  to  he  utilised  in  the  regions  of  small  relaxation  rates. 

The  time-marching  process  is  Initialised  by  defining  Atf(O),  ie  Aef(ci). 
Since  u(0)  must  he  zero,  it  follows  from  F.auation  13  that  Aef(ci)  is 
equal  to  df(0)/R.  This  quantity  is  readily  obtained  from  Reference  3. 

If  the  thermal  history  contains  temperature  discontinuities,  the  foregoing 
procedure  can  he  applied,  provided  that  the  strain  contributions  due  to 
the  discontinuities  are  dealt  with  separately. 

By  considering  the  net  change  in  the  circumferential  force  in  the  n  fibre 
layers  and  case  due  to  the  Vth  thermal  discontinuity  it  can  he  shown  that 
the  step  change  In  Aef(s),  at  C  “  Cki  Is  given  by 

,  „  xv  (ct  +  ofnA(0))(A(  AT(  %))) 

ACAe^))  -  ■  *  ’ 


A(0) 


Ef(0)tf 

Ectc 


(27) 


In  the  above  equation  the  quantities  A(Aef(g<))  and  ACATC^))  denote 
the  step  changes  in  Atf(c)  and  AT(^)  respectively  at  £  -  q<.  Once 
A(Aef(Ck))  has  i>ecri  calculated  the  value  is  added  to  the  immediately 
preceding  fibre  strain  change.  This  is  then  substituted  into  the  recurrence 
vector.  Equation  25,  and  the  solution  procedure  continued. 
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5  YIELD  INITIATION  PRESSL'FE 

The  composite  motor  mu_t  withstand  all  in-flight  loads  without  excessive 
distortion  or  failure.  Although  set-hack  forces  and  aeroheat  effects  may 
be  significant,  the  predominant  in-flight  load  will  always  be  the  firing 
pressure.  In  general,  motor  cases  have  been  designed  so  that  no  plastic 
deformation  occurs  under  this  pressure  loading.  Although  this  design 
approach  does  not  make  full  use  of  the  strength  reserves  available  in  the 
case,  and  in  this  sense  is  inefficient,  in  this  work  the  pressure  for  the 
onset  of  yield  has  been  adopted  as  the  basis  for  design. 

In  the  preceding  sections  the  stress  and  strain  components  for  the  fibres 
have  been  related  through  a  time-dependent  constitutive  eauation,  Equation 
4,  in  which  the  relaxation  modulus  is  known  to  decrease  with  time.  For 
loading  histories  extending  over  many  decades  (ie  10"^  h  t  10$  ;  the 
changes  in  the  relaxation  modulus  will  have  a  pronounced  effect  on  the 
stress  and  strain  states  in  the  motor.  For  very  short  duration  loading, 
eg  the  firing  pressure,  which  lasts  typically  for  less  than  one  second  and 
rarely  exceeds  several  minutes,  relaxation  effects  are  usually  neglected. 
Thus,  for  the  pressurisation  analysis,  the  stresses  and  strains  are 
related  through  the  conventional  Hookean  equations.  In  a  previous  analysis 
expressions  for  the  elastic  stress  changes  in  the  fibres  and  case 
corresponding  to  an  applied  internal  pressure  were  derived.  If  the 
pressurisation  process  is  regarded  as  an  instantaneous  load  application, 
these  expressions  can  be  readily  incorporated  with  Equations  18  and  19  to 
give  the  following  expressions  for  the  stresses  in  the  fibres  and  case 
at  t p,  the  time  of  pressure  appl ication: 

°efp(1,n)<tp>  -  Ef(o>|s9f(i,I,)(°>  +  i£f(?p)} 


^MiE-£2|c9f(i*n)(0) 


(28) 


+  Ae  (c')jdc' 

+  Pi  KA<0)  {i  -  [i-!i2]) 

t  f  ( 1+nA (0) )  ‘  2  1  P?  11 


,  i  «  1,2,... ,n. 


Oecp^V  ’  Ecf£0c(n)(°)  +  Aef(Cp) 
+  (of-ac)AT(?p)  } 


(29) 


+  piB  d  +  nvcA(0)  r1-Re2,l 
t  c(  1+nA  (0) )  1  2  1  ~R"ZJ 
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In  these  equations  vc  denotes  Poisson's  ratio  of  the  case  material,  Pp  is 
the  effective  throat  insert  radius3  and  pj  the  applied  internal  pressure. 

When  the  final  circumferential  stress  (Equation  29)  and,  where  appropriate, 
the  axial  stress  are  known,  the  pressure  for  the  onset  of  yield  can  be 
evaluated  using  the  von  Mises  yield  criterion,  the  most  appropriate  for 
the  materials  under  consideration'3. 

If  the  radial  stress  is  again  omitted  from  the  analysis,  the  mathematical 
form  of  this  criterion  is  given  by 

(oecp(")(cp))2  +  <ozcp<">)2 

"  0ecp(n)(V°zcp(n)  "  °y2 

(30) 

where  Oy  is  the  uniaxial  yield  stress  of  the  case  material,  and  ozct/n\ 
the  axial  stress  component  induced  in  the  case  during  firing,  is  given  by 


a  (n) 
°zcp 


(R2  -  Re2) 


Substituting  for  the  circumferential  and  axial  stress  components  in 
Equation  30  from  Equations  29  and  31,  a  quadratic  equation  for  the  pressure 
for  the  onset  of  yield  in  the  case  can  be  derived,  in  which  the  coefficients 
are  functions  of  the  winding  conditions,  fibre-relaxation  characteristics, 
end  conditions,  temperature  and  dimensions.  The  maximum  allowable  firing 
pressure  is  the  positive  root  of  this  quadratic. 


6  NUMERICAL  RESUITS 

To  illustrate  the  effects  of  fibre  relaxation  on  the  behaviour  of  the 
overwound  motor  case,  a  sample  motor  (the  same  as  that  considered  in 
Reference  3)  was  analysed  first  with  a  constant  uniform  temperature  equal 
to  the  winding  temperature  and  then  for  two  examples  of  varying  temperature 
conditions.  Motor  details  are  given  in  Table  1. 


Constant  Temperature  Conditions 


Because  the  reference  temperature  may  differ  from  the  winding  temperature 
the  following  analysis  is  presented  in  terms  of  the  reduced  time,  given 
by  5  ■  t/ax  (T),  rather  than  the  real  time  t. 
o 

For  the  purposes  of  illustration  the  viscoelastic  behaviour  of  the  overwind 
will  be  characterised  by  means  of  the  Dirichlet-Prony  series11  truncated 
after  the  first  two  terms. 


Eg ( C)  -  Ef(«)  +  (Ef(0)  - 


(32) 
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where  Ef(0)  and  F.f(®)  denote  the  Instantaneous  and  long-term  modulus 
values  respectively,  and  r  is  the  characteristic  time,  le  the  time  for  a 
modulus  change  of  633!  of  the  difference  between  the  Instantaneous  and 
long-term  values.  This  material  (the  standard  linear  solid)  displays 
both  instantaneous  and  asymptotic  long-term  elastic  behaviour,  and  is 
one  of  the  simplest  characterisations  of  a  viscoelastic  solid. 

With  such  a  material  it  is  possible  to  solve  Fquation  21  analytically  to 
obtain  the  time-dependent  fibre  strain  Aej(c).  Noting  that  AT(c)  is 
zero.  Equation  21  modifies  to: 


G(0)Aef(c)  ~  /  i£i£l£2Aef(c')dC' 

f  o  f 


(n) 

ec  (o)E.t_  [■ 


EfU) 


E{(0) 


1) 


(33) 


Substituting  from  Equation  32  into  Fquation  22  and  thence  into  Equation 
33,  a  first-order  linear  differential  equation  in  Atf(c)  can  be  derived, 
which,  with  the  boundary  conditions  A£f(0)  •  0,  has  the  solution  (see 
Annex  A). 


A ef(0  “  Aef(“)  {l-e'Y^t } 


(34) 


where  A  ef(»)  ■  eQc^T1^(0)(  1-1  )/(l+nA(  ®)  ) 

A(“)  *  Ef(«)tf/Ectc  , 
and  X  -  Ef (®)/Ef (0) 

Y  -  (l+nA(®))/(l+nA(0)) 


(35) 

(36) 

(37) 


(38) 

Once  Aef(c)  and,  as  a  consequence,  (see  Equation  16)  have  been 

obtained,  the  stresses  in  the  viscoelastic  fibres  and  the  case  resulting 
from  an  applied  internal  pressure  can  be  established.  For  the  fibres  the 
time-independent  and  time-dependent  stress  components  are  evaluated  by 
substituting  Equations  32  and  34  into  Equation  28.  Adding  the  resultant 
expression  to  the  pressure  component  (see  Equation  28)  gives,  after  some 
manipulation 
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°8fp(i,n)(cp>  "  feef(1,n>(0)  +  ^f(-)}Ef(!^) 

-  ief(-){Ef(0)e'YCP/T  + 

[Ef(0)-Ef  (»)]  [e“Cp/ /t]} 

1-1  P 


PjRACO) 

+  tf (I+nA(0)) 


,n. 


(39) 


Similarly,  for  the  case,  the  final  stress  is  evaluated  by  substituting 
from  Equation  35  for  Aef(“»)  in  Equation  34  and  thence  into  Equation  29 
to  give 


’Sep 


(n)(c„)  -  Fceec(n>(0)  „  +  (l-l)(l-e  YCP/'t) 


{1 


1+nA  (  “>) 


■} 


PiR 


tc(l+nA(0)) 


nvcA(0)  RP 
- - - F1--T 


]} 


(40) 


For  a  non-zero  axial  stress.  Equation  31  may  be  used  to  eliminate  p^  from 
Equation  40  in  favour  of  ozcp^n\  resulting  in  the  expression 


J8cp  -  Ec  eec<n>(0)(l  +  (A-l)(,l-el Vl)} 

1+nA ( ») 


(n), 


+  a  (n)  I  1  l  2 

°*co  > 


*cp  l(  l+nA(0)  ) 


+  vcnA(0)  j 


R  2 

(1-  !i_) 

R2 


(41) 


Once  the  total  circumferential  stress  in  the  case  is  known  the  pressure 
for  the  onset  of  yield  can  be  obtained  by  substituting  Equation  40  and, 
where  applicable.  Equation  31,  into  Equation  30.  Previous  work 3  has 
shown  that  the  quantity  Re  has  a  significant  effect  on  the  yield  pressure 
of  the  overwound  motor.  This  effect  is  examined  in  the  following  sections 
in  terms  of  the  maximum  and  minimum  values  of  Re. 

6.1.1  Open-Ended  Solution 

When  Re  approaches  R,  the  axial  stress  (Equation  31)  tends  to  zero  and 
the  circumferential  stress  in  the  case  (Equation  40)  becomes 
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’6cp 


(n) 


(Cp)  -  Eccec<n>(0){l 


+  (X-l)(l-e 


■V 


l+nA(«) 


+  PiR 

tc( l+nA(0) ) 


(62) 


Furthermore,  the  yield  Initiation  pressure.  Ply,  becomes  solely  dependent 
on  the  magnitude  of  the  circumferential  stress.  Hence,  from  Equation  62, 
it  follows  that 


piy 


(l+nA(0) )t- 


t°y_Ece0c 


(n)(0)[l+  }  • 

1+nA ( ») 


(63) 


For  the  Table  1  motor  the  variations  in  yield  pressure  with  n  for  several 
values  of  the  non-dimensional  reduced  time,  Cp  (iecp/t),  are  depicted 
in  Figure  2.  It  can  be  seen  from  the  figure  that  the  increase  of  p^v 
with  n  is  almost  linear  and  that  for  Cp“6,  p^y  becomes  almost  identical 
to  the  fully  relaxed  solution  (£p  ■  ®).  Since  n^(°)  (see  Reference 
2),  is  directly  proportional  to  the  winding  strain  t^pjy  varies  linearly 
with  Cw 


6.1.2  Closed-Ended  Solution 


When  the  quantity  Re  becomes  very  small,  ie  the  motor  approaches  its 
closed-ended  configuration,  the  axial  stress  (Equation  31)  becomes  non¬ 
zero.  By  substituting  the  two  principal  stress  expressions  (Equations  31 
and  60)  into  the  yield  criterion  (Equation  30)  a  general  expression  for 
the  yield  pressure  is  obtained.  For  the  Table  1  motor  the  variations  of 
the  yield  pressure  with  n  for  several  values  of  cp  are  shown  in  Figure 
3.  In  contrast  to  the  open-ended  solution  (Fig  2)  the  yield  pressure 
reaches  a  maximum  value  and  then  decreases  as  additional  layers  are 
Introduced.  For  the  instantaneous  elastic  solution  (ie  Cp  “^0),  this 
maximum  occurs  with  21  layers.  For  the  long-term  solution  ( Cp  “  “).  29 
layers  are  required.  It  can  be  shown  (see  Annex  B)  that  the  maximum 
yield  pressure  is  independent  of  the  initial  winding  configuration  and 
the  fibre-relaxation  characteristics,  and  is  equal  to  6oytc/RA.  For  the 
sample  motor  this  maximum  pressure  value  is  66  MPa. 

From  Figure  3  it  is  evident  that,  for  the  Table  1  motor,  the  yield  pressure 
decreases  with  increasing  Cp  for  all  n  values  less  than  26  and  Increases 
with  increasing  Cp  for  all  n  values  greater  than  26.  These  contrasting 
trends  are  further  evident  in  Figure  6,  where  the  yield  pressure  is 
plotted  against  cp  for  different  values  of  X,  for  the  two  arbitrary  n 
values  of  12  and  60.  For  n  ■  12  the  yield  pressure  is  seen  to  decrease 
with  decreasing  X,  whereas  for  n  »  40  the  opposite  occurs.  It  is  noteworthy 
that  when  Ef(~)  becomes  zero  the  residual  stresses  in  the  case  disappear, 
and  the  pressurisation  response  of  the  motor  will  depend  only  on  the 
instantaneous  elastic  modulus  of  the  fibres. 
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An  Insight  Into  these  effects  can  he  obtained  from  Figure  5,  in  which  the 
elliptical  von  Mises  yield  boundary  for  a  bi-axial  stress  state  is  plotted; 
the  line  bB  in  the  figure  represents  the  instantaneous  elastic  solution 
of  Equation  41.  It  has  been  shown  previously  3  for  elastic  conditions 
that,  when  the  prestrained  motor  is  pressurised,  the  stress  response  proceeds 
along  the  line  bB  until  the  point  V  is  reached,  at  which  point  yield 
occurs.  The  corresponding  yield  pressure  is  proportional  to  the  clcp(nVc^ 
ordinate  at  this  point  and  can  be  evaluated  using  Equation  31.  More 
generally,  for  viscoelastic  conditions,  the  q^  ^n'(0)/ay  value 
representing  the  length  Ob  is  given  by  the  first  term  on  the  right-hand 
side  of  Equation  41,  which  is  time-dependent,  and  the  slope  of  the  line 
is  obtained  from  the  second  term  of  Equation  41,  which  is  time-independent. 
With  relaxation,  therefore,  the  stress  response  line  bB  moves  to  the 
right  in  the  figure,  retaining  the  original  slope,  with  the  changing 
length  Ob*  given  by  Ogcp  n' ( Cp)/ <3y  •  The  yield  intercept,  V,  becomes 
V*  as  indicated.  It  is  clear  that,  depending  upon  the  initial  position 
of  V  in  relation  to  the  acme  of  the  ellipse  and  the  maximum  shift  of  the 
bB  line,  the  pressure  for  first  yield  may  increase,  decrease  or  increase 
to  a  maximum  and  then  decrease  with  (p. 

6.2  Varying  Temperature  Conditions 

In  the  previous  section  constant  temperature  conditions  were  assumed  and 
it  was  possible  to  obtain  stress  solutions  and  yield  pressure  values  in 
closed  form. 

By  contrast,  for  varying  temperature  conditions  the  stresses  are  obtained 
numerically  by  solving  the  Integral  Equation  21  for  Lcf(0  and  substituting 
into  Equation  29.  In  representing  the  stress  behaviour  graphically, 

Equation  31  is  used  to  eliminate  p^  from  Equation  29  to  give  a  6imple 

linear  relationship  between  Ogcp^^Cp)  ant*  dfecp^  oE  the  form 

‘0cp(n)<«p>  "  Ecle9C(n)<0)+ief(«p)+(«1f-^>iffS)} 


+  o  <">{  1 

*cp  l(l+nA(0)) 


+v  nA(D))  j  J 


(44) 


For  the  purpose  of  illustration,  the  above  equation  can  be  represented 
exactly  as  in  Figure  5.  The  slope  of  the  bB  line  is  obtained  from  the 
second  term  on  the  right-hand  side  of  Equation  44,  and  is  independent  of 
both  temperature  and  time.  The  intercept  Ob*  is  given  by  the  first  term 
on  the  right-hand  side  of  Equation  44;  it  will  he  dependent  on  both 
temperature  and  time  and  its  variation  will  be  governed  by  the  relative 
magnitude  of  the  two  strain  components  and  (af-ofc  )il(  Cp). 
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To  illustrate  the  effects  of  varying  temperature  conditions  upon  the 
yield  pressure,  the  Table  1  motor  has  been  studied,  overwound  with  12  and 
40  fibre  layers  and  subjected  to  temperture  histories  of  the  form  AT(t)  - 
a  sin  ui t,  first  with  a  ■  50°C  and  id  »  2 it  (hours- *)  and  then  with  a  • 

50°C  and  id  “  0.2s  (hours-1).  The  motor  was  assumed  to  be  closed-ended 
with  Re  equal  to  zero.  Five  complete  temperature  cycles  were  considered 
for  both  histories.  Three  forms  of  material  behaviour  for  the  overwind 
have  been  assumed; 

a.  elastic, 

b.  viscoelastic  with  temperature-independent  properties 

(see  Section  2.1). 

c.  viscoelastic  with  "thermorheologically  simple"  properties 

(see  Section  2.2). 

Nolte  and  Findley  16  found  that  a  generally  varying  temperature  history 
could  be  satisfactorily  represented  by  a  series  of  time  intervals  over 
which  either  the  temperature  or  the  temperature  rate  was  assumed  constant. 

In  the  present  work,  however,  only  the  latter  will  be  considered.  For 
both  temperature  histories  each  thermal  cycle  was  approximated  as 
illustrated  in  Figure  6. 

6.2.1  12-layer  Configuration 

For  the  12-layer  configuration  the  intersection  of  the  bB  line  with  the 
ellipse  (point  V,  Fig  5)  lies  to  the  right  of  the  acme  of  the  von  Mises 
yield  ellipse.  Consequently  the  yield  pressure  will  decrease  with 
decreasing  case  pre-stress  (as  indicated  by  the  length  Ob),  and  vice  versa. 
Since  of  is  less  than  a,.,  it  follows  for  the  elastic  conditions  that 
the  yield  pressure  increases  with  increasing  temperature,  the  opposite 
occurring  for  a  decreasing  temperature.  For  the  Table  1  motor  this 
behaviour  is  illustrated  in  Figures  7  and  8,  where  the  yield  pressure  is 
plotted  against  time  for  the  two  temperature  histories  considered. 

When  viscoelastic  effects  are  Included  these  trends  are  modified.  The 
yield  pressure  may,  depending  on  the  temperature  history  and  fibre- 
relaxation  characteristics,  increase  or  decrease  with  increasing  temperature. 
The  following  features  are  evident  in  Figures  7  and  8: 

a.  In  Figure  7,  with  the  faster  temperature  cycle,  the 
variations  in  pjy  derived  from  the  temperature-independent  analysis 
(curve  2)  follow  the  temperature  changes  with  a  reduction  in 
successive  peak  values  due  to  relaxation.  It  is  evident  that  the 
time  scales  of  the  temperature  changes  and  the  relaxation  effects 
are  such  that  initially  the  former  dominates.  Relaxation  effects 
are  seen  to  be  practically  negligible  after  5  cycles. 

b.  By  contrast,  for  the  thermorheologically  simple  analysis 
(curve  3,  Fig  7),  it  is  seen  that  after  an  initial  Increase,  Piy 
decreases  with  Increasing  temperature,  indicating  that  relaxation 
effects  quickly  become  dominant.  The  effect  is  reversed  before 
the  first  temperature  peak,  Implying  near- completion  of  the 
relaxation  process. 
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c.  For  the  slower  temperature  cycle  (Fig  8),  it  is  clear 
that  relaxation  effects  are  dominant  over  the  first  quarter-cycle 
for  the  temperature-independent  analysis  (curve  2),  since  pjy 
decreases  with  Increasing  temperature*  (This  is  in  contrast  with 
curve  2,  Figure  7  -  see  a.  above.)  There  are  no  further  relaxation 
effects  after  the  first  half-cycle. 

d.  Apart  from  the  absence  of  the  initial  increase  in  pjy, 
there  is  little  difference  in  the  prediction  derived  from  the 
thermorheo  gically  simple  analysis  for  the  fast  and  slow 
temperature  cycles  (see  curve  3,  Fig  8  and  curve  3,  Fig  7). 

6.2.2  40-layer  Configuration 

In  contrast  with  the  12-layer  configuration,  the  bB  line  (Fig  5)  for  the 
40-layer  configuration  intersects  the  von  Mises  yield  locus  to  the  left 
of  the  acme  of  the  ellipse  and  the  yield  pressure  will  increase  with 
decreasing  ogcp^n ^(Cp)/o„  and  vice  versa.  It  follows  therefore  that 
the  trends  described  in  Section  6.2.1  will  be  reversed.  This  is  evident 
in  Figures  9  and  10.  Additional  features  worth  highlighting  are: 

a.  There  is  a  marked  lack  of  symmetry  of  the  variations  in 
Ply  about  its  mean  value,  particularly  for  the  elastic  analysis, 
for  both  long  and  short  period  temperature  variations. 

b.  The  amplitude  of  the  p^y  variations,  again  particularly 
for  the  elastic  analysis,  is  considerably  greater  than  that 
observed  in  Figures  7  and  8. 

c.  In  Figure  10  for  the  slower  temperature  cycle,  the 
predicted  piy  values  from  the  two  viscoelastic  analyses  are 
practically  identical  after  the  first  temperature  cycle. 

7  CONCLUSIONS 

1  Solution  techniques  and  illustrative  results  have  been  presented 
for  the  stresses  developed  in  a  cylindrical  metallic  motor  case 
circumferentially  overwound  with  a  prestrained  viscoelastic  fibre  for 
both  constant  and  varying  temperature  histories.  The  internal  pressure 
required  for  the  initiation  of  yield  in  the  reinforced  case  has  also  been 
determined. 

2  Relaxation  effects  may  result  in  either  an  increase  or  a  decrease 
in  the  yield  initiation  pressure. 

3  With  a  constant  temperature  history  and  assuming  open-ended 
conditions  In  the  cylinder  (le  motors  with  an  effective  nozzle  radius 
which  approaches  that  of  the  motor  case)  the  yield  initiation  pressure 
for  a  particular  example  Increases  almost  linearly  with  the  number  of 
fibre  layers  (n)  and  decreases  non-linearly  with  time  for  a  given  n, 
approaching  the  fully  relaxed  conditions  at  times  greater  than  about  four 
times  the  characteristic  time  of  the  fibre- 
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4  With  a  constant-temperature  history  and  closed-ended  conditions 
In  the  case,  the  yleld-lnltlatlon  pressure  increases  to  a  maximum  with 
Increasing  n  and  thereafter  decreases  as  additional  layers  are  Introduced. 
There  Is  a  markedly  non-linear  time-dependence,  which  varies  with  n. 

5  The  variations  In  yleld-lnltlatlon  pressure  resulting  from  two 
sinusoidally  varying  temperature  histories  have  been  studied  for  the 
closed-ended  configuration.  The  effects  on  the  yleld-lnltlatlon  pressure 
of : 


a.  the  period  of  the  temperature  cycle 

b.  the  material  behaviour  assumptions 

for  two  distinct  winding  configurations  have  been  illustrated  and  discussed. 
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9  NOTATION 

aT  (T) 
o 

A(0) 

A(«») 

E(t) 

Ec 

EfU) 

EfJ 

f(c) 

G(?) 


n 

Pi 

R,Re 

SJ.L 

t 

lf  >lc 


shift  function  (see  Equation  7) 

see  Equation  27 

see  Equation  36 

relaxation  modulus  at  time  t 

Young's  modulus  of  case  material 

relaxation  modulus  of  fibre  overwind  at  reduced 
time 

jth  Prony  coefficient  (see  Equation  23) 
see  Equation  22 
see  Equation  22 

number  of  Prony  coefficients  associated  with 
Equation  23 

number  of  applied  fibre  layers 
internal  pressure 

mean  and  effective  throat  insert  radius  of 
case  respectively 

recurrence  vector  (see  Eauation  25) 
real  time 

effective  fibre  thickness  and  case  thickness 
respectively 
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AT(?) 

u(c) 

Of.Oc 

Y 

6f(c),6cU) 

e(t ) 

e0f(1’n)(c).  Eqc^H?) 
AEf  (c),  Aec(c) 

C 

A 

C 

X 

o(t) 

o0f(1'n>Cc),o9c(n>(;) 

Oy 

o  <"> 

°zcp 

T 

Tj 

MC 

(i) 

Subscripts 

P 

TOT 

AT 


shift  function  reference  temperature 

time-dependent  uniform  temperature  change 

time-dependent  radial  expansion  of  motor  case 

thermal  expansion  coefficients  of  fibre  over¬ 
wind  and  case  respectively 

see  Equation  38 

time-dependent  interference  displacement  of 
fibre  and  case  respectively  (see  Equations  10 
and  11) 

time-dependent  strain  (see  Equation  1) 

time-dependent  circumferential  rtrain  components 
in  the  ith  fibre  layer  and  case  respectively 

time-dependent  strain  changes  in  fibre  and  case 
(see  Equations  13  and  If) 

reduced  time  (see  Equation  6) 

non-dimensional  reduced  time  (ie  q/ t) 

fibre  modulus  ratio  (ie  Ef ( »)/Ef (0) ) 

time-dependent  circumferential  stress  (see 
Equation  1) 

time-dependent  circumferential  stress  components 
in  the  it*'  fibre  layer  and  case  respectively 

case  yield  stress 

axial  stress  component  in  metallic  case 

characteristic  time  of  standard  linear  solid 

relaxation  time  associated  with  j1*1  Prony 
coefficient  (see  Equation  23) 

Poisson '8  ratio  of  metallic  case 

angular  frequency 


due  to  pressure  application 
total 

associated  with  a  uniform  change  in  temperature 
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TAB  IF  1  Sample  Composite  Motor  Case  Data 


Case  Young's  modulus: 

Ec 

SK 

72.4  GPa 

Fibre  relaxation  moduli: 

Ef(0) 

* 

93.1  GPa 

Ef  (“) 

46.5  GPa 

Characteristic  time  of  fihre: 

1 

- 

l.Ohr 

Poisson's  ratio: 

^c 

C 

0.3 

Case  yield  stress: 

°y 

- 

470.0  MPa 

Wall  thicknesses: 

me 

2 .  Onm 

tf 

X 

0. 1mm 

Motor  radius: 

p 

X 

- 1 

47 . 0mm 

Winding  strain: 

K 

4.0  x  lO”3  1 

Thermal  expansion  coefficients: 

<v 

K 

22.5  x  10-&/-C 

°f 

X 

3.6  x  UT*/°C 
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FIG. 
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FIG.  5 
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FIG.  7  P-  vs  t_  FOR  aT(»  =  50"C  sin  2wt( n  =  12  CLOSEP-ENDED) 
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FIG.  8 
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FIG.  8  P-  vs  t  FOR  AT(t)  =  50  C  sin  0  27it  (n  =  1 2  CLOSEP-ENDED) 
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FIG.  9 


o 

in 


o 


o 

CO 


CO 

.c 

a 


o 

c\i 


O 
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FIG.  9  P;„  vs  t„  FOR  &T(t)  =  50  C  sin  2u\  (n  =  40  CLOSEP-F.NPEP) 
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ANNEX  A 

ANALYTICAL  SOLUTION  TO  THE  CONSTANT  TEMPERATURE  ANALYSIS  INTEGRAL  EQUATION 
(EQUATION  T5J 

In  Section  6.1  the  Integral  equation 

G(0)Aef(5)  "  /'  C' 


E9c(n)(0)Ectc[ 


Ef  ( 0 
Ef(0) 


-1} 


was  derived  for  the  constant  temperature  analysis  where,  from 
Equation  22, 

G(?)  ”  nEf ( c)tf  +  Ectc 


(Al) 


(A2) 


For  the  standard  linear  solid  the  relaxation  modulus  Ef({)  appearing 
in  Equations  Al  and  A2  can  be  readily  characterised  by  the  truncated 
Dirichlet-Prony  series^,  in  the  form 

Ef(c)  -  Ef(-)  +  Eie-5/T 

(A3) 


in  which  (see  Equation  32) 
EX  -  Ef(0){l-X} 


(A4) 


With  this  characterisation  Equation  A2  modifies  to: 
G(C)  -  ntf{Ef(-)+E1e-c/T}  +  Ectc 


(A5) 


Substituting  from  Equations  A4  and  A3  for  Ef({)  in  Equation  Al  and 
noting  from  Equation  A5  that 


<*GU-t')  -  G,e~<?-i;')/T 
dC’  1 


where 


Gl  -i 


(A6) 


(A7) 


Al 

UNLIMITED 


IWIIMJTFD 


the  following  modified  integral  equation  is  obtained 

C(0)Aef(O  -  Cj'e^C-C’^Aeffc’JdC' 
o 


■0C 


<n) 


(°)Fctc 


,F.f  («)+Eie 

1  moT 


-Jj  T 


-1} 


(AB) 


Differentiating  the  above  equation  with  respect  to  C  gives,  after  some 
manipulation,  the  following  first-order  linear  differential  equation  in 
Aef  ( c) 


<*Aef(c)  +  ,1 

dC  't 


£l 

G(0) 


}itf  ( 0 


£9c(n>(0)(  X-l) 
t(  l+nA(  0) ) 


(A9) 


which,  with  the  boundary  condition  Aef(0)-C,  has  the  solution 
Aef(c)  -  Aef(-)(l-e-Yt/T) 


(A10) 


with 


Acf(-) 


c  (n)(o)_QiLL 

^  1+nA  ( ■») 


(All) 


and 


Y  a  l-HiA(-) 

l+nA(0)  (A12) 

Once  Acf(c)  has  been  obtained  the  corresponding  strain  change  for  the 
case  material  Aec(0  then  follows  from  Equation  16.  The  time-dependent 
fibre  and  case  strains  and  stresses  are  then  evaluated  via  Equations  12 
and  14  and  Equations  IB  and  19  respectively. 
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ANNEX  B 


MAXIMUM  YIELD  PRES  SUP  E  ANALYSIS 
ISOTHERMAL  CONDITIONS 


If  the  radial  stress  component,  which  Is  a  second-order  quantity. 
Is  omitted  from  the  von  Mises  yield  criterion.  Equation  30,  the 
simplified  expression  is  obtained 


<°9cp<n)>2  +  “ 


zcp 


30cp 


“zcp 


(n) 


(B1 ) 


where  the  circumferential  and  axial  stress  components,  and 

qzcp<n>  respectively,  have  the  functional  form 

°ecP(n)  '  °ecp<n){Piy^} 


“zcp 


(n) 


“zcp 


(n) 


{Piy> 


(B3) 


These  quantities  are  also  functions  of  the  winding  conditions,  moduli, 
geometry,  etc,  but  for  the  purposes  of  this  work  they  are  assumed  to  take 
prescribed  values* 


Differentiating  Equations  Bl,  B2  and  B3  with  respect  to  tp  and  setting 
3piy/3Cp  to  zero,  the  following  relationships  for  the  maximum  yield 
pressure  conditions  can  be  derived 


either 


or 


3oecp 

Kp 

?o0cp 


(n) 

(n)  - 


0  • 


°zcp 


(n) 


0 


(BA) 

(B5) 


Since  O0Cn  Is  a  monotonically  increasing  function  of  £p  (see 
Equation  40)  it  follows  that  only  the  second  condition.  Equation  BS, 
applies.  Equation  B5  can  now  he  used  to  eliminate  o^p^11  from 
Equation  Bl  to  give  the  simple  relationship 


o  (n){p 
zcp  iymax 


iL_  o 

/3  X 


(BS) 


where  Piymax  denotes  the  maximum  yield  pressure.  Using  Equation  B6, 
the  maximum  yield  pressure  then  follows  from  a  simple  rearrangement 
of  Equation  31,  ie 


Bl 

UNLIMITED 


UNLIMITED 


4t  £  OyR 

piymax  "  /3(R2_Rg2) 


(B7) 


Once  Piymax  has  Bean  obtained,  the  time  to  the  maximum  yield  pressure 
condition,  ?pmax»  f°r  the  standard  linear  solid  then  follows  by  substituting 
from  Equation  B5  for  OgCpCn)  in  Equation  41,  and  thence  from  Equation 
B6  for  oZCp^n)  in  the  resulting  equation  to  give,  after  some  manipulation 


max 


-Iln 


Y 


Oy(24rl) 

/3EcAcf  (») 


+ 


e8c(n)(0)  +1  i 
1  1 


where 


Atf  (”) 


eac^^OKl-l) 

l+nA(«0 


,  (see  Equation  35), 


(B8) 


(B9) 


TT+zKmr  r 


+  v  nA(0)  ] 
c 


Cl-  ) 


(BIO) 


Equation  B8  can  now  be  easily  solved  for  5pmax"  Only  positive  values 
are  admissible. 
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